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Dark Matter as Dense Color Superconductor 
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Canada 

We discuss a novel cold dark matter candidate which is formed from the ordinary quarks during the QCD 
phase transition when the axion domain wall undergoes an unchecked collapse due to the tension in the wall. If 
a large number of quarks is trapped inside the bulk of a closed axion domain wall, the collapse stops due to the 
internal Fermi pressure. In this case the system in the bulk, may reach the critical density when it undergoes a 
phase transition to a color superconducting phase with the ground state being the quark condensate, similar to 
BCS theory. If this happens, the new state of matter representing the diquark condensate with a large baryon 
number B > 10 20 becomes a stable soliton-like configuration. Consequently, it may serve as a novel cold dark 
matter candidate. 
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1. Introduction 

The presence of large amounts of non-luminous 
components in the Universe has been known for a 
long time. In spite of the recent advances in the 
field |Q, the mystery of the dark matter/energy 
remains: we still do not know what it is. The 
main goal of this talk is to argue that the dark 
matter could be nothing but well-known quarks 
which however are not in the "normal" hadronic 
phase, but rather in some "exotic", the so-called 
color superconducting (CS) phase. 

This is a novel phase in QCD when light quarks 
form the condensate in diquark channels, and it 
is analogous to Cooper pairs of electrons in ordi- 
nary superconductors described by BCS theory. 
There existence of CS phase in QCD represents 
our first crucial element for our scenario to work. 
The study of CS phase received a lot of atten- 
tion last few years S. It has been known that 
this regime may be realized in nature in neutron 
stars interiors and in the violent events associ- 
ated with collapse of massive stars or collisions of 
neutron stars, so it is important for astrophysics. 
The goal of this talk is to argue that such con- 
ditions may occur in early universe during the 
QCD phase transition, so it might be important 
for cosmology as well. 

The force which squeezes quarks in neutron 
stars is gravity; the force which does a similar job 
in early universe during the QCD phase transition 



is a violent collapse of a bubble formed from the 
axion domain wall. If number of quarks trapped 
inside the bulk is sufficiently large, the collapse 
stops due to the internal Fermi pressure. In this 
case the system in the bulk may reach the critical 
density when it undergoes a phase transition to 
CS phase with the ground state being the diquark 
condensate. We shall call the configuration with a 
large number of quarks in color superconducting 
phase formed during the QCD phase transition 
as the QCD-ball. Therefore, an existence of the 
axion domain wall represents our second crucial 
element for our scenario to work. As is known, 
the axion is one of the favorite candidates for the 
cold dark matter, see recent reviews ||. In the 
present scenario the axion field plays the role of 
squeezer rather than dark matter itself. 

As will be shown below, once such a configu- 
ration is formed, it will be extremely stable soli- 
ton like particle. The source of the stability of 
the QCD-balls is related to the fact that its mass 
Mb growth as Mb ~ B 8 / 9 for the large baryon 
(quark) charge B and becomes smaller than the 
mass of a collection of free separated nucleons 
with the same baryon charge. The region of the 
absolute stability of the QCD-balls is determined 
by inequality mjv > d Q B B which can be always 
satisfied for very large B > B c . However, even 
for sufficiently large B but smaller than B c , i.e. 
1 < B < B c , we still expect that QCD-balls 
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will be very long-lived metastable states. This is 
because the elementary excitations carrying the 
baryon charge in CS phase and hadronic phase 
are very different: in CS phase they carry color 
charge along with baryon charge; in hadronic 
phase they are color singlets. Therefore, if suf- 
ficiently large number of quarks is trapped in- 
side the axion bubble during its shrinking, it may 
result in formation of an absolutely stable ( or 
long lived, metastable) QCD-ball with the ground 
state being a diquark condensate. Such QCD- 
balls, therefore, may serve as the cold dark matter 
candidate and contribute to Qdm — 0.3. 

Strictly speaking, the QCD-balls being the 
baryonic configurations, would behave like non- 
baryonic dark matter. In particular, QCD-balls, 
in spite of their QCD origin, would not contribute 
to fish 2 ~ 0.02 in nucleosynthesis calculations 
because the QCD-balls would complete the for- 
mation by the time when temperature reaches the 
relevant for nucleosynthesis region T ~ XMeV. 
Once QCD-balls are formed, their baryon charge 
is accumulated in form of the diquark condensate, 
rather than in form of free baryons, and in such 
a form the baryon charge is not available for nu- 
cleosynthesis. 

2. QCD-balls 

Crucial for our scenario is the existence of a 
squeezer, axion domain wall which will be formed 
during the QCD phase transition. We assume 
that the standard problem of the domain wall 
dominance is resolved in some way as discussed 
previously in the literature ||. We also assume 
that quarks which are trapped inside of the ax- 
ion bubble, can not easily escape the bulk when 
the bubble is shrinking. In different words, the 
axion domain wall is not transparent due to the 
QCD sandwich structure of the wall as discussed 
m [@ i H ■ The collapse is halted due to the Fermi 
pressure. Therefore, we assume that a large num- 
ber of quarks remains in the bulk of volume V 
when the Fermi pressure cancels the surface ten- 
sion. Exactly this condition determines a typical 
radius Ro of the QCD-ball at the moment of for- 
mation. The pressure due to the surface tension is 
P a = while the fermi pressure Pf in the bulk 



can be easily estimated from the thermodynam- 
ical potential for the non-interacting relativistic 
fermi gas, Q = —PfV with Q given by: 
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Here fj, is the Fermi momentum of the system to 
be expressed in terms of the fixed baryon charge 
B trapped in the bubble; 
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g ~ 2N c Nf = 18 is the degeneracy factor for 
massless degrees of freedom. The condition P a ~ 
Pf determines a typical radius of the QCD-ball 
with a fixed initial baryon charge B at the mo- 
ment of formation: 
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Here a ~ fam^f^ is the axion domain wall ten- 
sion with f a ~ (10 10 - 10 12 )GeV being con- 
strained by the axion search experiments. As the 
first approximation, we neglected many contribu- 
tions in ([lj ||) which might be important numeri- 
cally but can not drastically change our main re- 
sults. In particular, we neglected the quark-quark 
interaction on the Fermi surface, which brings the 
system into superconducting phase for relatively 
large \i > AOOMeV B. The corresponding con- 
tribution gives some correction to the energy of 
order (^-) 2 with A ~ IQQMeV is the supercon- 
ducting gap, and can be neglected for the asymp- 
totically large fj,. We also neglected the so-called 
bag constant contribution which is a phenomeno- 
logical way to simulate the confinement. Para- 
metrically this contribution is also suppressed as 
( Aq ^ D ) 4 k 11 ^ in reality might be important. Let 
us emphasize: we are not attempting to solve the 
problem of equilibrium between CS and hadronic 
phases^]. Rather, we are attempting to estimate 

1 Such a problem has been recently discussed in rcf. 
where the interface region between nuclear matter an. 
CFL (color flavor locking) superconducting phase was an 
alyzed. 
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the initial size Rq of the QCD ball which is formed 
due to the collapse of the axion domain wall 
halted by the fermi pressure. Other contribu- 
tions, such as quark-quark interaction or the bag 
constant, may change the numerical estimates for 
Rq however they can not replace the main player 
of the game, the fermi pressure which is respon- 
sible for the stopping of the collapse. 

Once the QCD ball with radius i?o is formed, 
the bulk of this object being in CS phase will be 
quite stable with respect to the decay on the free 
nucleons if some conditions are met, see below. 
The stability remains intact even if the original 
axion domain wall would decay: the configura- 
tion remains stable due to the CS phase itself, 
and not because of the support from the axion do- 
main wall. The only role the axion domain wall 
plays in our scenario is the squeezing the mat- 
ter at the very first instant in order to bring the 
large number of quarks into the CS phase. The 
later stages of the evolution of the QCD ball will 
determine the precise structure of this configura- 
tion at the equilibrium, including the interface re- 
gion separating the CS phase with nuclear matter 
phase; we do not expect, though, that this evolu- 
tion would considerably change our estimates for 
Rq presented above. 

Once a typical radius i?n for a QCD ball is es- 
timated (||), the total energy E ( fermion energy 
+ axion domain wall energy) and baryon number 
density n of the QCD-ball can be also estimated 
with the result 

E ~ 4wi? 2 + gV^ ~ l^c^B^a 1 ^ 



The most important consequence of this result is 
the behavior E ~ B 8 / 9 which implies that for suf- 
ficiently large B the QCD-ball becomes an abso- 
lutely stable object. It happens when > 
i.efl 

B> B c ~ (2 9 7rc 2 -^ r ) 3 ~ 10 33 , n c ~ 9n , (4) 

m N 

where no is the nuclear saturation density, no ~ 
0.16(/m)~ 3 . Numerically, B c ~ 10 33 corresponds 

2 Our normalization for the baryon charge corresponds to 
B = 1 for the quark, thus factor ^p-- 



to stabilization radius R c ~ W 12 GeV 1 and mass 
E c ~ 10 33 GeV. Here we use /„ ~ 10 w GeV which 
corresponds to a ~ 1.8 • 10 8 GeV 3 . As we men- 
tioned in Introduction, we expect that the QCD- 
balls with much smaller B can live long enough 
to serve as a dark matter. Therefore, in what fol- 
lows we use a more realistic value for B^ xp ~ 10 20 
which follows from the experimental bounds on 
masses and fluxes of nontopological solitons, see 
below. 

Few remarks are in order. First, we note that 
there is an upper limit on value B above which 
the QCD-balls can not be formed at the very 
first instant during the collapse of the domain 
walls. This follows from the expression (Q) for 
the baryon number density which decreases when 
B increase. Therefore, at some point, the density 
n reaches the magnitude where the phase transi- 
tion between nuclear matter and CS phase occurs. 
Numerically it happens at B ~ 10 2 B C . However, 
it does not preclude from formation of a larger 
object at the later stages of the QCD-ball evolu- 
tion. 

As our next remark, we would like to mention 
some similarity between the QCD-balls (which is 
the subject of this letter) and Q-balls[j7| which 
are nontopological solitons associated with some 
conserved global Q charge. In both soli- 
ton mass as function of Q has behavior, similar 
to our eq. (^), and therefore, it may become a 
stable configuration for relatively large Q charge. 
Therefore, an effective scalar field theory with 
some specific constraint on potential (when Q ball 
solution exists) is realized for QCD in high den- 
sity regime. The big difference, of course, is that 
underlying theory for QCD-balls is well known: it 
is QCD with no free parameters, in huge contrast 
with the theory of Q-balls. Formal similarity be- 
comes even more striking if one takes into account 
that the ground state of the CFL phase in QCD 
is determined by the diquark condensate with the 
following time dependence (\I r \f r )* ~ (e* 2/lt ), with 
if? being the original QCD quark fields, and /x be- 
ing the chemical potential of the system, see for- 
mula (40) from rcf. |||. As is known, such time- 
dependent phase is the starting point in construc- 
tion of the Q balls 0. 

Finally, we want to constraint the QCD-ball pa- 
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rameters from available data. We follow || and 
assume that a typical cross section of a neutral 
QCD-ball with matter is determined by their ge- 
ometrical size, ttRq. We also assume that the 
QCD-balls is the main contributor toward the 
dark matter in the Galaxy, in which case their 
flux F should satisfy 



F < 



Pdmv 



AttM 



B 



7.2-10 5 %^cm- 2 sec- 1 sr- 1 , (5) 
Mb 



where Pdm 1S the energy density of the dark 
matter in the Galaxy, p DM ~ 0, ^f , and v ~ 
3 • 10~ 3 c is the Virial velocity of°the QCD-ball. 
We identify Mb in the expression (||) with the 
total energy E (||) of the QCD ball at rest with 
given baryon charge B. Different experiments 
discussed in |J give the following bound on the 
flux of neutral soliton-like objects: F < const. ■ 
10 _16 cm _2 sec _1 sr _1 , This translates to the fol- 
lowing lower limit of the neutral QCD-ball mass 
Mb and baryon charge B, 

M e B xp > 2 • 10 21 GeV, B exp > 1.6 • 10 20 . (6) 

These experimental bounds are well below the 
critical line of the absolute stability of the QCD- 
balls, given by eq. (||). However, as we mentioned 
above, we expect that QCD-balls will be very 
long-lived metastable states even if their baryon 
charge is below the critical limit (^). Therefore, 
we consider experimental constraint (^|), rather 
than a theoretical stability bound (Q) as a more 
realistic bound on Mb, B for the QCD-balls. 

We have seen that the observed relation ~ 
£Idm within an order of magnitude finds its nat- 
ural explanation in this scenario: both contri- 
butions are originated from the same physics at 
the same instant during the QCD phase transi- 
tion. As is known, this fact is extremely difficult 
to explain in models that invoke a dark matter 
candidate not related to baryons. We have also 
seen that without adjusting any parameters the 
baryon density inside the QCD-ball falls exactly 
into appropriate region of the QCD phase dia- 
gram where color superconductivity takes place. 
The scenario is no doubt lead to important con- 
sequences for cosmology and astrophysics, which 
are not explored yet. In particular, some unex- 
plained events, such as Centauro events, or even 



the Tunguska-like events (when no fragments or 
chemical traces have ever been recovered) , can be 
related to the very dense QCD balls. If this is 
the case, the arrival directions should correlate 
with the dark matter distribution and show the 
halo asymmetry. Also: recent discovery []Io| sug- 
gests that the matter in some stars could be even 
denser than nuclear matter. It could be also re- 
lated to the very dense QCD balls. Therefore, the 
"exotic", dense color superconducting phase in 
QCD, might be much more common state of mat- 
ter in the Universe than the "normal" hadronic 
phase we know. In conclusion, qualitative as our 
arguments are, they suggest that the dark matter 
could be originated at the QCD scale. 
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